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CN ; Abstract 

We explore the possibility to observe hard exclusive three-jet production in early 
^ I LHC runs, corresponding to diffractive dissociation of the incident proton into three 

jets with large but compensating transverse momenta. This process is sensitive to 
^ ■ the proton unintegrated gluon distribution at small x and to the distribution of the 

three valence quarks in the proton at small transverse distances. The corresponding 



Dh' 



^! cross section is calculated using an approach based on kt factorization. According to 

r-| I our estimates, observation of hard diffractive three-jet production at LHC is feasible 

for jet transverse momenta g_L ~ 5 GeV. 

>' 

cn ■ !• The physics potential of forward detectors at LHC, within and beyond the standard 

model, is attracting a lot of attention, cf . [H [21 El HI E] . In this Letter we explore the possibility 
^ i to observe hard exclusive diffractive dissociation of a proton into three hard jets in proton- 

(X) ! proton collisions 

O ■ p{p^) + p(p2) -^ jet{qi) + jet{q2) + jet(g3) + p{p'^) , (1) 



X 



cf. [6] In this process one proton stays intact and the other one dissociates into a system of 
three hard jets separated by a large rapidity gap from the recoil proton, see Fig. 1. The main 
c^ \ aim of our study is to estimate the cross section of this reaction and the corresponding event 

rates at LHC and Tevatron. 

Note that we are interested in exclusive three-jet production which constitutes a small 
fraction of the inclusive single diffraction cross section. The exclusive and inclusive mechanisms 
have different final state topologies and can be distinguished experimentally. A characteristic 
quantity is e.g. the ratio Rjets of the three-jet mass to the total invariant mass of the system 
produced in the diffractive interaction. Exclusive production corresponds to the region where 
Rjets is close to unity. This strategy was used recently at the Tevatron [7] where central 
exclusive dijet production, pp — *• p + jet+iet + p, in double-Pomeron coUisions was measured 
for the first time. 

Exclusive dijet production in the central region has much in common with the exclusive 
Higgs boson production process, pp -^ p + H + p. In [8J it was argued that studies of exclusive 
dijet production and other diffractive processes at the early data runs of the LHC can provide 



Figure 1: Proton dissociation into three jets. The unintegrated gluon distribution includes the hard 
gluon exchange as indicated by the dashed square. 

valuable checks of the different components of the formalism. Indeed, this was the main moti- 
vation for Tevatron experiment. The exclusive 3-jets production in single diffraction (Q offers 
another interesting example since factorization of hard and soft interactions in this case is less 
complicated. In particular, the fluctuation of a proton projectile into a state with small trans- 
verse size, which is the underlining mechanism for ([1]), suppresses secondary soft interactions 
that may fill the rapidity gap. Thus one can get an access to the gluon distribution at small 
a; in a cleaner environment, having no problems with gap survival probability and factoriza- 
tion breaking that introduce major conceptual theoretical uncertainties in the calculations of 
diffract ive Higgs production. 

Our approach to exclusive three-jet production derives from experience with coherent pion 
diffraction dissociation into a pair of jets with large transverse momenta which was measured 
by the E791 collaboration [HI HO]. The qualitative features of the E791 data have confirmed 
some earlier theoretical predictions [HI [121 [13] • ^ strong A-dependence which is a signature 
for color transparency, and a ~ 1/q^ dependence on the jet transverse momentum. These 
features suggest that the relevant transverse size of the pion r± remains small, of the order of 
the inverse transverse momenta of the jets r±_ ~ l/q±. 

On a more quantitative level, we have shown [HI [THj that coUinear factorization is violated 
in dijet production due to pinching of singularities between soft gluon (and quark) interactions 
in the initial and final state. However, the nonfactorizable contribution is suppressed com- 
pared to the leading contribution by a logarithm of energy so that in the double logarithmic 
approximation Ing^ Ins/g^ collinear factorization is restored. Moreover, to this accuracy hard 
gluon exchange can be "hidden" in the unintegrated gluon distribution JF(x, q±). Thus, in the 
true diffraction limit, for very large energies, hard exclusive dijet production can be considered 
as a probe of the hard component of the pomeron. The same interpretation was suggested 
earlier in [16] within the kt factorization framework (see also [H]). The double logarithmic 
approximation turns out to be insufficient for the energy range of the E791 experiment, but 
might be adequate for the LHC In this Letter we present an estimate for the cross section for 
the reaction ([T]) based on the generalization of these ideas. 

2. At leading order the jets are formed by the three valence quarks of the proton, see 
Fig. 1. We require that all three jets have large transverse momenta which requires at least 



two hard gluon exchanges. One of them can be effectively included in the high-momentum 
component of the unintegrated gluon density (the bottom blob) as indicated schematically by 
the dashed square, but the second one has to be added explicitly since the hard pomeron only 
couples to two of the three quarks of the protonl^ The necessity for an additional hard gluon 
exchange makes calculation of proton diffraction dissociation more difficult as compared to the 
meson case. 

Our notation for the momenta is explained in Fig. [H We neglect power corrections in 
transverse momenta of the jets and also proton and jet masses so that Pi = P2 = p'i = Q'l = 
oi — Qz — 0- The jet momenta are decomposed in terms of momenta of the initial particles 

qk = OikPi + PkP2 + qk±, k = 1,2,3 (2) 

where 

qi± + q2± + q3± = 0, ai + 02 + "3 = 1 , Pk = q^x/i^ks) , (3) 

The three-jet invariant mass is given by 

M' = {q^+q, + q,f = ^ + ^ + ^, ( = — = (3, + (3, + (3, . (4) 

«! 0:2 as s 

where s = {pi+P2)^ = 2pi-p2 is the invariant energy. Assuming that the relevant jet transverse 
momenta are of the order of 5 GeV, the typical values of the ( variable at LHC are in the 
range C ~ 10"^ ^ 10"^ 

At high energies, an amplitude is predominantly given by its discontinuity in the s-channel 
which usually implies that the amplitude is almost purely imaginary. In our case the situation is 
more complicated since in the physical region of ([T]) the amplitude develops a cut in the variable 
M^ as well, and it remains complex even after taking the s-channel discontinuity. Nevertheless, 
the s-channel discontinuity of the amplitude can be expressed, in our approximation, in terms 
of the unintegrated gluon distribution JF(x, k±). 

The relevant Feynman diagrams in leading order of perturbative QCD are shown in Fig. 2. 
They can be divided into three groups which differ by the attachments of the t-channel gluons 
to the quark lines (shown by crosses). In diagrams (f)-(j) the hard gluon exchange takes 
place between quark q2 and q^. Hence the transverse momentum of one of the t-channel 
gluons coincides with that of the quark gi. As a consequence, the contribution of this group of 
diagrams involves the unintegrated gluon distribution at the same scale, J-'{C: qi±)- ^^ diagrams 
(k)-(o) the hard gluon exchange connects the quarks gi and q^, the transverse momentum of the 
i-channel gluon coincides with the momentum of the quark ^2 and, therefore, the unintegrated 
gluon distribution enters at this scale, J-'(C, ?2±)- Similarly, the contribution of diagrams (a)-(e) 
is proportional to ^(C, q3±)- 



^Alternatively, one can consider double-pomeron exchange in the t-channel. This contribution is suppressed 
by a power of the jet transverse momentum so it is of higher twist. 
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Figure 2: The leading-order contributions to proton disintegration into three jets within kt 
factorization. The points where the t-channel gluons are attached to the quarks are shown by 
crosses. 



Accordingly, we have three different contributions to the amplitude: 
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where / Da' = J^ da'ida'2da'^6{l — Y^ a'j) corresponds to the integration over the quark momen- 
tum fractions in the incident proton, e^^^ describes the color state of the final quarks, N = 3 
is the number of colors. The dimensionless quantities Ci are expressed in terms of different 
Dirac structures where for convenience we introduce a "positron-like" Dirac spinor v, 



{u{q2)) =Cv{q2), {u{q2hi,) = -C-i^v{q2) . 



(6) 



Here C is the charge-conjugation matrix. The diagrams with the 3-gluon vertex do not con- 
tribute due to vanishing color factors. Therefore we need to calculate in total 12 nontrivial 
diagrams (4 for each group). 

The calculation is straightforward, though rather tedious. Here we present the final results 
only: 

Cf^j = [V u{qi)^iv{q2) u{q3)fi2l5N{pi) - ^M(gi)^i 75^(92) M(g3)^2A^(Pi)] 

ai{a3 + a'^) — ai(a3 + 0:3) 



X 



+ [Vu{qi):^2v{q2) u{q3)^2'lbN{pi) - Au{qi)^2'lbv{q-2) u{q3)^2N{pi)] 

, ^2/^2 - Cdpi aipi + a^ps + (^3 - a'2)p2 

X 7 ; 7: — ^^Tr — TTT^ — : — 7r'^ 7: — : — rr + 



a3"3/33S^][a2(/^2 + Ps) - ai(3i + ie] [-«2"3/32S^][a3(/32 + Ps) - "lA + ^e] 
"2/^2 + "3/^3 



Q;2Q;3/92/93(a;2 + as)s^ ^ 

+ [Vu{qi)^i-f5v{q2) u{q^)'^2N{pi) - Au{qi)^iv{q2) u{q2)^2'lbN {pi)] 

-ai{a2 + Oi'2) ai{a2 + OL. 

X 7 ; 7, ^TTv 7T7, : 7n 7, : rr + 



«2"3/52S^][tt3(/^2 + /^s) - ttiA + ^e] [-a3a3/53-s^]["2(/^2 + P:i) - "lA + ^e] 
+ \yu{qi)]^2lbv{q2)u{q3)^2N{pi) - Au{qi)]^2v{q2)u{q:i)'^2l^N{pi)] 

a'2J32 + ail3i -oiil3i + q;3/93 + ("3 - tt2)/^2 



+ 



-a3a3/33s3][a^(/32 + 13^) - ai(3i + ie] [-a2"3/92S^][a3(/52 + Ps) - "lA + ie] 
a'2p2 + a'^Ps 



a'2cy'^l32l33{a2 + as)s^ ^ 
+ T [u{qi)v{q2) u{q3)'y5N{pi) - u{qi)'^^v{q2) u{q3)N{pi)] 

. ^3 -'^1 ("2 - ai)/92 + («2 - «3)/^3 

X 7 m — ;tt7 — 7T7 — : — 7n 7\ — : — rr + 



a^/?2s2] [a'^ {(32 + Ps) - ai A + ie] [-a'^Pls^] [«^ {P2 + Ps) - aipi + ie] 
{a2 -ai){a2p2 + a3p3) 



a2a'^P2Pi{oi2 + as)^^ 
+ T [^(gi) {ipif,P2,,cr^^) v{q2) uiqshbNipi) - u{qi) {ipif,P2u<y t,u) l^v{q2) u{q-i)N{pi)] 



2(1 - a'2) ^ 2((ai + a2)p2 + (^2 - a'^)Pz) 



-a'2p2s'^] [a'^ {p2 + Ps) -aipi + ie] [-a'^Pls^] [a'2 {P2 + Ps) - aipi + ie] 
2{ai + 0:2) (0:2/32 + a'^Ps) 



a'2a'^P2Pl{,Oi2 + 03) s^ 



(7) 



La^e = [V u{qi):^iv{q2) u{q3):^2l5N{pi) - Au{qi)fii-f5v{q2)u{q3)fi2N{pi)] 
^ , 2K - «i) ^ 2K - "2) 



-a'2(32S^] [a[{Pi + P2) - aa/^s + ie] [-a'^pis^] [a'^{pi + P2) - as/^s + ie] 
+ [V u{qi):^2v{q2) u{q3)^2l5N{pi) - Au{qi)^2l5v{q2)u{q3):^2N{pi)] 
Wi , 2/^2 



+ 



-a'2/32s3] [a[{Pi + /32) - 03/33 + le] [-a\pis^] [a^(/3i + P2) - asPs + ^e] 



a[a2PiP2{cti + 02)3^ ^ 
+ T [u{qi)v{q2) u{q3)'y^N{pi) - u{qi)'y5v{q2) u{'l3)N{pi)] 

-2^3 , 2^3 



-a'2/?2s2] [ai(/5i + /52) - asf33 + ie] [-a[f3is^] [a'^if^i + (32) - a^fi^ + ie] 

The functions ^(a'j^, ctg, 0:3), V(a']^, 02, 03) and T(a']^, 0:2, 0:3) are the leading-twist hght-cone 
nucleon distribution amphtudes defined as in p^ 

(0|e*^V„(aiz)M^^(a2^)<(a3^)|iV(j9i)) = 
= V{^^C)^p (75iV(pi)), + A ihl^C)^^ (iV(pO), + T (tcr.^p'^C)^^ {Yl,Nip^)l^ , (9) 

where z"^ = and a^,^ = (z/2)[7^,7j. Each invariant amphtude V,A,T depends on the 
scalar products aiPi ■ z and can be represented as the Fourier transform of the corresponding 
distribution amplitude, e.g. 

n.....)^/B«'e-^i:-;"^y(„:,. (10) 

Finally, iZfc-o is given by the expression similar to Cf-j in Eq. ([7]) with the replacements 
a\ ^-> 0:2, a\ '^ a'2 and /3i <->• /32 everywhere except for the arguments of the distribution 
amplitudes. 

One remark is in order. Another possible mechanism for the exclusive proton disintegration 
into three jets could be the exchange of three hard gluons in the t- channel. Such a contribution 
could involve two different color structures, proportional to j°^'^ and d°^'^ ^ which correspond to 
different C-parity in the t— channel. We found by explicit calculation that the color factor 
~ yabc corresponding to the C-parity-even exchange vanishes. The C- parity odd contribution 
~ d'*^^ is related to odderon exchange and presumably small. 

3. The differential cross section can be written as 

I Alp dai(ia2(ia3(5(l -ai -^2 -tts) ,2- j2- j^ji, /hn 

^^ = ovo ^8 2 ^ ^1^ q2dtd(\)t 11 

2^(27r)*s^ a\a2az 

where t = {j>2~ p'2)^ is the Mandelstam t variable of the pp scattering and ^t is the azimuthal 
angle of the final state proton. In our kinematics, for large transverse momenta of the jets 
and small t, one can neglect effects of azimuthal correlations between the jets and the final 
proton. Hence d^t integration is trivial and gives a factor 27r. For the t dependence we assume 

6 



a simple exponential form, da/dt ~ e^*, and use 6 ~ 4 -^ 5 GeV^ for the slope parameter which 
is a typical value which describes HERA data for hard exclusive processes: DVCS and vector 
meson electroproduction at large Q^. Thus, the integration over the proton recoil variables 
gives a factor J dtdcpt — > ^• 

Since our calculation is only done to double logarithmic accuracy, we use the simplest model 
for the unintegrated gluon distribution as given by the logarithmic derivative of the usual gluon 
parton distribution xgi^x^Q"^) 

d 

•^(a;, gi) = 7T1 2^9{x,q]_). (12) 

omq]_ 

The numerical estimates presented below are obtained using the CTEQ6L leading-order gluon 
distribution as provided by [19j. We also used the simplest, asymptotic model for the nucleon 
distribution amplitude 

V{a[) = T{a[) = 120 fNa[a'^a'^ , ^(a^ = . (13) 

The normalization parameter /jv is scale-dependent. To leading-logarithmic accuracy 

/„(,).;„(,„) (1^)"^, ,14) 

where /3o = 11/3A^ — 2/3nf. The existing QCD sum rule estimates 

rjiMfNifx = 1 GeV) = (5.0 ± 0.3) x 10"^ GeV^ [20] , 

m^/^(;x = 1 GeV) = (5.1±0.4) X lO^^GeV^ [2l] (15) 

are somewhat larger compared with a very recent n/ = 2 (unquenched) lattice calculation 

/jv(/i = 2GeV) = (3.14±0.09) X lO^^GeV^ [22]. (16) 

The given number corresponds to the lattice spacing a ~ 0.067 fm; a continuum extrapolation 
was not attempted. For definiteness we use the value f^ = 5.0 x 10~^ GeV^ at 1 GeV as input, 
and evolve it to the relevant scale. 

The integration over the phase space of the three jets was done numerically, restricting the 
longitudinal momentum fractions to the regioro 

0.1 < ai,a2,a3 < 0.8 (17) 

and requiring that the transverse momentum of eac/i jet is larger than a given value go = Q'±,min- 
For the value qo = 5 GeV we obtain for the integrated three-jet cross section at the LHC energies 

^LHC _4 . f fNJqo) y f as{qo) Y f 5GeV y (1R\ 

^^-^•-^-^P^-Ur-lO-^GeV^J lT2rJ [-^) ■ ^^^) 



^We use this rather conservative cut condition to assure a clear three-jet event selection. 



Assuming the integrated luminosity for the first LHC runs in the range 100 pb^^ to 1 fb^^ 
an observation of this process at LHC seems to be feasible. Note that the effective power 
o" ~ l/^Q (fitted in the go = 3 ^ 8 GeV range) is somewhat stronger than the naive power 
counting prediction a ~ l/^o- This effect is due to the strong ^ dependence of the unintegrated 
gluon distribution: larger values of go imply larger invariant masses M ^ of the three-jet system 
(jl]) and consequently larger C, = M'^/s. The sizeable cross section for go = 5 GeV is in 
fact an implication of the expected rise of the LO gluon distribution more than two times 
as ( is decreasing by roughly a factor of 50 when going from Tevatron to LHC. The existing 
parameterizations of the LO gluon distribution at C ~ 10~^ differ from each other by ~ 30%. 
The unintegrated gluon distribution fll2p enters as a square in the prediction for the cross 
section, therefore, the study of exclusive three-jet events at LHC may provide a valuable 
constraint for the gluon distribution at small momentum fractions. 

A comparison of the three-jet exclusive production at LHC and the Tevatron can be es- 
pecially illuminating in this respect since other uncertainties do not have significant impact 
on the energy dependence. For Tevatron kinematics, assuming the value g_Lmin = 3 GeV, our 
estimate for the cross section (fitted in the range go = 2 -=- 4.5 GeV) is 

o^r - 5opb . ( .,{"'*> V (mx f— V. (19) 



^''' "^ V4.7-10-3GeV2>' V 0.255 



% 



Note that in this case M^ ~ 100 GeV^ and C ~ 10 ^ ^ 10 ^ where the gluon distribution 
is much better known: The typical difference between existing parameterizations is of order 
~ 10%. 

In the most naive approximation, the longitudinal momentum fraction distribution of the 
jets is expected to follow that of the valence quarks in the proton: The momentum fraction 
distribution of a jet, arbitrary chosen in each event, is proportional to the proton distribution 
amplitude squared 

^r^ [ Da'6{a~a[)\M(^')\\ (20) 

da I 



where 0Ar = V — A |T8]. In reality, a hard gluon exchange leads to a certain redistribution of 
the longitudinal momenta so that the resulting a-dependence is more complex, see e.g. [HI [15]. 
To illustrate this effect, in Fig. [3]we compare our calculated normalized jet momentum fraction 
distribution (averaged over quark flavors) to the dependence in ( 120]) : The two distributions are 
similar, but the one resulting from the QCD calculation is shifted towards lower momentum 
fractions compared to the simple dependence in Eq. (!20l) . Note that measurement of the 
valence quark momentum fraction distribution in a pion presented the main motivation for the 
E791 experiment [9l [10] which, in turn, triggered detailed studies of such reactions in pQCD. 

4. To summarize, in this Letter we have studied the exclusive diffractive dissociation of a 
proton into three jets with large transverse momenta in the double-logarithmic approximation 
of perturbative QCD. This process is interesting in the broader context of diffractive processes 
at LHC, which will be studied using forward detectors, and in particular can be used to 
constrain the gluon distribution at very small values of Bjorken x. According to our estimates. 
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Figure 3: The normalized longitudinal momentum fraction distribution of the jet (solid curve). 
For comparison, the dependence corresponding to Eq. (1201) is shown by dashes. 

an observation of such processes in the early runs at LHC is feasible for jet transverse momenta 
of the order of 5 GeV. 
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